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We present a new dimensional scaling transformation of the Schrodinger equation for tiie two elec- 
tron bond. Tills yields, for the first time, a good description of the two electron bond via D-scaling. 
There also emerges, in the large-D limit, an intuitively appealing semiclassical picture, akin to a 
molecular model proposed by Niels Bohr in 1913. In this limit, the electrons are confined to specific 
orbits in the scaled space, yet the uncertainty principle is maintained because the scaling leaves in- 
variant the position-momentum commutator. A first-order perturbation correction, proportional to 
1/D, substantially improves the agreement with the exact ground state potential energy curve. The 
present treatment is very simple mathematically, yet provides a strikingly accurate description of the 
potential energy curves for the lowest singlet, triplet and excited states of H2. We find the modified 
D-scaling method also gives good results for other molecules. It can be combined advantageously 
with Hartree-Fock and other conventional methods. 



Quantum chemistry has achieved excellent agreement 
between theory and experiment, even for large molecules, 
by using computational power to overcome the difBculty 
of treating electron-electron interactions Here 
we present a new version of an unconventional method 
to treat electronic structure |E IE Sill- This emulates an 
approach developed in quantum chromo dynamics , by 
generalizing the Schrodinger equation to D dimensions 
and rescaling coordinates Q- 

Early work found the tutorial D-scaling procedure of 
Witten Pi can be dramatically improved; the ground 
state energy of He was obtained accurate to 5 significant 
figures by interpolation between the D ~ 1 and D —^ 00 
limits ^ , and to 9 figures by a perturbation expansion in 
1/D 10]. However, the scaling procedure which worked 
well for atoms [51161 did not prove successful for two- 
center problems 0, 0; e.g., for H2 that procedure did 
not yield a bound ground state, (see our Fig. 

In our present approach, the large-D limit makes con- 
tact with the Bohr model of the H2 molecule In 
this way we obtain, for the first time, a link between pre- 
and post-quantum mechanical descriptions of the chemi- 
cal bond (Bohr-Sommerfeld vs Heisenberg-Schrodinger). 
Marked improvement is achieved by including the lead- 
ing correction term in l/D and a rudimentary adjust- 
ment of the D-scaling. Fig. ^ shows potential energy 
curves for H2 obtained with our simple approach. Dots 
comprise a synthesis of experimental data and computa- 
tions employing many terms in variational wavefunctions 
[1^. Our simple method gives surprisingly accurate re- 
sults and holds promise for numerous applications. 

We first outline our method as apphed to Ho and 
then indicate how it differs from what preceded 0, Q. 
Fig. 13 displays electron distances in the H2 molecule. 
All distances are expressed in terms of the Bohr length 
Go = /me^, where m is the electron mass, and energies 
are in the Hartree unit e^/oo. We start with H'^ = E'^, 
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FIG. 1: Potential energy (solid curves) of the ground and a 
few excited states of H2 obtained from the Bohr model with 
D-scaling analysis. Dots are the "exact" energies [l^ . The 
inserted figures on the right hand side depict the two nuclei 
of charge Z and Bohr's "planetary" orbits for the electrons in 
states (see also Fig. 0J. Dashed curves are 
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FIG. 2: Electronic distances in H2 molecule. The nuclei A 
and B are fixed a distance R apart. 



The Coulomb potential energy V is given by 

Z Z Z Z 1 
+ — + ^, 1 

Val ni ra2 n2 ^12 R 

in terms of distances defined in Fig. |21 In cylindrical 
coordinates 
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where R is the internuclear spacing and (f> the dihedral 
angle between the planes containing the electrons (i = I, 
2) and the internuclear axis. 

We proceed by endowing each vector with D carte- 
sian coordinates 0. The potential energy V is retained 
|Q| in the three dimension form of Eq. whereas the 
Laplacians in the kinetic energy take the form 
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We then scale coordinates by and energy by 1//^, 
with f = {D — l)/2, and transform the wavcfunction "if 

by 



This recasts the Schrodinger equation as 
{Ki+ K2 + U + V)<S> = E<i>, 

where 
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In the limit D — > 00 the derivative terms in Ki are 
quenched. The corresponding energy Eoc for any given 



internuclear distance R is then obtained simply as the 
extremum of the effective potential, U + V , given by 

E = I- ( ^ + ^] +V{pi,p2,zi,Z2,<j),R). (6) 
2 \Pi P2J 

This is exactly the energy function that applies to the 
Bohr model of the molecule [lll |. 

The usual D-scaling procedure 0, 1^ involves set- 
ting up the full Laplacian in D-dimension and trans- 
forming the wavefunction by incorporating the square 
root of the Jacobian via 'J J~^^^^, where J = 
{piP2)^~^{sin(j>)^~^. Then, on scaling the coordinates 
by /2 and the energy by I//2, the Schrodinger equation 
in the limit D 00 yields 
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+ V{pi,p2,Zi,Z2,(l),R), (7) 



which differs from Eq. |(HJ by the factor l/sin2 0. 

Our procedure, designed to reduce to the Bohr model 
at the large- 13 limit, instead incorporates only the radial 
portion of the Jacobian in transforming the wavefunc- 
tion via Eq. This has important consequences. Fig. 
13 displays the D 00 potential energy curve of Eq. 
JT)) (dashed curve, "full- J"); which exhibits no binding. 
However, our "Bohr model" limit obtained from Eq. © 
yields a good zero-order approximation for the ground 
state (curve 2 in Fig. OJ. It is surprisingly accurate at 
both large and small internuclear distances R. Also, the 
model predicts the ground state is bound with an equi- 
librium separation i?e = 8/(9 — V3) ~ 1-10 and gives 
the binding energy as Eb ~ 3(2 — \/3)/8 a.u.« 0.100 
a.u.= 2.73 eV. The Heitler-London calculation (shown in 
Fig. ^ dashed curve), obtained from a two-term vari- 
ational fimction, gives Re = 1.51 and Eb = 3.14 eV 
(22], whereas the "exact" results are Re — 1.401 and 
Eb = 4.745 eV Q. 

For the triplet ^Y,^ state, as seen in Fig. ^ the Bohr 
model energy function of Eq. ® gives a remarkably close 
agreement with the "exact" potential curve and is in fact 
much better than the Heitler-London result (which, e.g., 
is 30% high at i? = 2). 

In essence, D-scaling procedures resemble gauge trans- 
formations. Many varieties of scaling are feasible, sub- 
ject only to the constraint that as D ^ 3 the scaled 
Schrodinger equation reduces to the correct form. The 
basic aim is to devise a scaling that removes the ma- 
jor, generic D-dependence, enabling the easily evaluated 
D 00 limit to approximate the D = 3 energy. With the 
"full- J" scaling previously used 0, when D is increased 
the (sin0)-°~'^ factor in the Jacobian forces cj) towards 
90°, while minimization of electron-electron repulsion re- 
quires (f> 180°. The effect is to overweight electron 
repulsion; this is the chief source of the failure to obtain 
chemical bonding in previous work. Our new procedure 
avoids such overweighting by retaining the D — 3 form 
for the (/)-part of both the Jacobian and the Laplacian 
of Eq. |2Jl. Thereby (f) remains a fully quantum variable 
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FIG. 3: Energy E{R) of H2 molecule in the limit 73 — > 00 
calculated from Eq. JTJ (dashed curve) and from the Bohr 
model of Eq. Q (solid curves). Curve 1 corresponds to a 
symmetric configuration obtained by Bohr Tl^ and pictured 
in Fig. |31(top). Curve 2 describes an asymmetric solution (not 
found by Bohr), see Fig. |1| (bottom). Lower solid curve is the 
improved ground state E{R) after including 1/D correction. 



as -D — > 00, rather then being converted to a semiclassi- 
cal parameter along with the p and z coordinates. This 
much improves description of the electron repulsion and 
hence the chemical bonding. 

The scaling procedure enables, in the largc-D limit, cal- 
culations to be carried out in the scaled space that are 
entirely classical. The extremum equations dE/dz = 
and dE/dp = are equivalent to Newton's second law 
applied to the motion of each electron. Respectively, 
they specify that the net Coulomb force on the electron 
along the z— axis vanishes and that the projection of the 
Coulomb force perpendicular to the molecular axis bal- 
ances the centrifugal force. Although the electrons are 
thereby confined to specific orbits in the scaled space, 
the uncertainty principle is nonetheless satisfied. This is 
so because the conjugate momenta are scaled inversely to 
the coordinates, leaving the position-momentum commu- 
tator invariant. The continuous transition between the 
scaled space and the unsealed space in effect relates clas- 
sical trajectories at large-D to corresponding quantum 
distributions at D=3. This aspect becomesparticularly 
evident when treating electronic tunneling Q- 

Fig. ^ displays the "exact" electron charge density 
along the molecular axis in the ground state of H2 for 
internuclear spacing R — 0.8 and 1.4 a.u. Circles show 
electron orbits in Bohr's model. The orbit positions for 
any R actually coincide with the maxima in the charge 
density. This provides a link between the wave mechani- 
cal and Bohr (D — > 00 limit) treatments of the H2 bond. 

The ground state E{R) can be substantially improved 
by use of a perturbation expansion in powers of 1/-D, de- 
veloped by expanding the effective potential of Eq. ® 




FIG. 4: Distribution of the electron charge density in the H2 
molecule along the molecular axis z. The nuclei are fixed a 
distance R apart. Circles are electron orbits in Bohr's model. 



in powers of the displacement from the minimum |3 ; for 
He this has yielded highly accurate results Terms 
quadratic in the displacement describe harmonic oscil- 
lations about the minimum and give a l/D correction 
to the energy. A symmetry breaking point occurs at 
Rc — 1.2, beyond which the electron orbits move apart 
(c.f. Fig. 0J. Such symmetry breaking is a typical fea- 
ture exhibited as Z or R is varied at large- Z3 jT^ ,1^] . The 
1/D correction works well at points substantially below 
or above Rc- Results for those regions thus can be com- 
bined. This involves transforming the axial coordinates 
to zi ± Z2, in order to separate the double- well structure 
that occurs in zi — Z2- With the other coordinates fixed at 
their values at the minimum oiU + V , a one-dimensional 
Schrodinger equation is solved to take into account the 
double- well mode. This contribution to the 1/D correc- 
tion corresponds to electron resonance or exchange. The 
result gives good agreement with the "exact" E{R) over 
the full range of R (lower solid line in Fig. |3Jl. The 
1/D correction predicts the equilibrium separation to be 
Re — 1.38 with binding energy Eb = 4.50 eV. 

The Bohr and D-scaling techniques taken together hold 
promise for numerous applications. In particular, these 
provide a new approach to treating excited states. For 
example, in our analysis the energy of ls2s state of the He 
atom is obtained as an extremum of the energy function 
E = n\/2r\~\-n\/2r\ + V{Y'i_,r-2), where ni — 1 and 712 = 
2; ri, r2 are electron radius vectors and V is the Coulomb 
potential energy. This yields the value of —2.159 a.u. 
which differs by 0.7% from the "exact" ls2s energy of 
—2.144 a.u. For other excited states of He as well as 
more complex atoms the combination of the Bohr and 
D-scaling approaches also provides accurate results; we 
will discuss this elsewhere. 

Fig. ^ demonstrates application of our technique to a 
few excited states of the H2 molecule. In treating 
excited states, we incorporate D-scaling analysis at large 
R and the exact E{R) of the H^ molecular ion which 
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FIG. 5: Ground state E{R) of HeH, He2 and BeH molecules 
calculated within the Bohr model (solid curves). The HeH 
curve is shifted down for clarity. 

provides a good description in the remaining region. We 
have also found the present — oo limit (Bohr model) 
gives good results for other molecules; examples so far 
treated include HeH, He2, and BeH, pictured in Fi g. El 
and LiH, Li2, Be2, and the triatomics BeH2 and H3 [Tq . 



Another useful strategy is to combine the present ap- 
proach with conventional electronic structure methods. 
At D — 3, evaluation of the correlation energy, i^corr 
(error in the Hartree-Fock approximation) is the major 
difficulty. However, at D ^ 00, Econ can be evaluated 
exactly. Results for He and other atoms show that 
-Ecorr for D — > c» is smaller than but comparable to 
that for D = 3. For the ground state of H2 we find an 
accurate energy curve E{R) can be obtained by adding 
the D — > 00 correlation energy to the E{R) given by the 
Heitler-London effective charge method. The result is 
practically identical to the curve obtained from the 1/D 
correction (Figs. ^andlSl). 

Our modified D-scaling procedure reincarnates the 
Bohr model. This requires only elementary concepts 
and (laptop) computations yet provides a rather good 
description of electron-electron interaction and chemi- 
cal bonding. The procedure is readily applicable to 
many-electron molecules, both ground and excited states. 
These results encourage efforts to further improve D- 
scaling and to augment conventional variational methods 
for electronic structure to incorporate the exact correla- 
tion energy attainable at the large-D limit. 
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